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1 INTRODUCTION 



SU(2) lattice gauge theory with dynamical fermion at non-zero chemical potential 
and at finite temperature is studied. We focus on the influence of chemical potential 
for quark condensate and mass of pseudoscalar meson at finite temperature. Hybrid 
^J • Monte Carlo simulations with Nf = 8 staggered fermions are carried out on 12 X 

^y. ' 12 X 24 X 4 lattice. At /3 = 1.1 and rriq =0.05,0.07,0.1, we calculate the quark 

condensate and masses of pseudoscalar meson consisting of light and heavier quarks 
for chemical potential fi = 0.0,0.02,0.05,0.1,0.2. 

o 

^j. • The study of strong interactions at finite baryon density lias a long history 

f^ I Recent theoretical speculation is that the ground state of dense baryonic mat- 

f^ ■ ter may. he more exotic El' 13. Nambu-Jona- Lasinio type and instanton liquid 

^P I models □' B suggest that the QCD vacuum at sufficiently high baryon density 

could become a color superconductorlj. 

Considerable progress has been made in simulating lattice quantum chro- 
modynamics ( LQCD ) at nonzero temperature Q. But, there has been lit- 
(j^ I tie progress in understanding the theory at nonzero chemical potential even 

rS^ • though a pioneering work was presented sixteen years agoQ. The basic difficulty 

in simulation of LQCD at finite density is that the effective action becomes 
. . - complex due to the introduction of chemical potential. Standard algorithrns 

rN I such as hybrid molecular dynamics ( HMD )l3 or hybrid Monte Carlo ( HMC )lI 

j^ ' can not be applied in this situation. In order to avoid the complex determinant, 

early simulations Jjsed quenched approximation, but they do not correspond 
to the real world a. In this situation, a possible way is to examine response 
of physical quantities to chemical potential oLm, = 0. Quark- number suscepti- 
bilities have been studied by S.Gottlieb et alu and give us the signals for the 
chiral- symmetry-restoration phase transition. Study of response of hadron 
masses to chemical potential is now under investigation by QCD — TARO 
collaborationllj. 

Another way is to study LQCD for color SU(2) system where action re- 
mains real. Hamiltonian approach at zero temperature and finite density has 
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been studied by E.B.Gregory et aall, and gives us some results for vacuum en- 
ergy, chiral condensate, baryon-number density and its susceptibility. Numer- 
ical simulation at non zero chemical potential has been performed by S. Hands 
et aO, by using several algorithms; HMD algorithm ; HMC algorithm ; Two- 
Step Multi-Boson ( TSMB ) algorithmtj on rather small lattices ( mostly on 
6"* lattice ) at zero temperature. They studied relation of chiral condensate, 
pion mass and chemical potential at /x > ^ as well as /i < '^. 

In this paper, we focus on the effects of chemical potential to masses and 
condensate at finite temperature ( /x=0.02, 0.05, 0.1 and 0.2 ) on larger lattices 
12'^ X 4 and 12^ x 24 x 4. One of interests is precise study of masses of meson with 
nonzero chemical potential since QCD Sum Rule suggests linear dependence 
on density. Ej Firstly, we shall observe chiral condensate with nonzero chemical 
potential near critical temperature. Secondly, we shall examine meson masses 
consisting of light-light quarks and light-heavy quarks. The organization of the 
paper is as follows: In Section 2 we state our procedure. Simulation is described 
in section 3. We show some preliminary results of the present simulation in 
Section 4. Summary and conclusions are stated in Section 5. 

2 STAGGERED FERMION ACTION AT NONZERO CHEMI- 
CAL POTENTIAL 

At non zero chemical potential, staggered fermion matrix M is given by 

3 

+ r?4(x)[e''C/.,45,,,+4 - e-^C/J,4(5^,,_4], ^^^ 

where U is SU(2) link variables and fj, is chemical potential, ry^ is staggered 
phase and m is quark mass parameter. 

The matrix M satisfies the following relation 

The path integral for the partition function is given by 

Z = J[DU]det{M)^f/^e^s= , .„. 

where Sg is plaquette action written as , 

Sy^PY.[^^\Tr{Up + Ul)], 

p (4) 



where B — ^r . 

In cases of simulations by HMC algorithm, we need positive definiteness of 
hamiltonian. Since det(M) is real for SU(2) system, the positive definiteness 
condition is achieved by taking det(M)^ which includes n/ = 8 dynamical 
fermions. Note that the standard even-odd partitioning of the fermion matrix 
which is commonly used to reduce the number of flavors can not be applied for 
simulations at finite chemical potential. Thus n/ = 8 is the minimum number 
of flavors for SU{2) HMC simulations with dynamical fermions. The present 
study is performed with Uf = 8. 

3 SIMULATIONS 

Firstly, we have done simulation on 6^ lattice and confirmed that our results are 
consistent with the previous work by S. Hands et alc3. After that, numerical 
simulations on 12^ x 4 and 12^ x 24 x 4 lattices have been carried out for 
parameters summarized in Table 1. 



Tabic 1: List of the data sample 



Np 


TTlq 


Ai 


lattice size 


/3 


8 


0.05 
0.07 
0.1 


0.0 , 0.02 

0.05 , 0.1 

0.2 


6x6x6x6 


0.1-2.5 ( for test ) 


12 X 12 X 12 X 4 


1.0-2.0( to determine /3c ) 


12 X 12 X 24 X 4 


1.1 ( for simulations ) 



Number of configuration for every parameter is also given in Table 2. 







Table 2: Statistics of the data 








Ai = 0.0 


fi = 0.02 


// = 0.005 


Ai = 0.1 


Ai = 0.2 


ruq 


= 0.05 


120 


180 


160 


80 


30 


niq 


= 0.07 


140 


170 


160 


160 


30 


niq 


= 0.10 


180 


170 


160 


140 


70 



Most of the simulations are carried out on small chemical potential and 
at fixed /3 = 1.1 which is little below (3c- Time step of HMC algorithm is 
dt=0.01. After 1000 trajectories of thermalization, data are taken at every 100 
trajectories. Since our simulation is restricted in a region of small chemical 
potential, the algorithm works without difficulty. But, necessary CPU time 
increases with the chemical potential. A typical trend is shown in table 3. At 
II = 0.2, average CPU time is 3.5 times longer than that at /^ = 0.0. 



Tabic 3: CPU time of getting one configuration at /3 = f .1, rriq = 0.05 
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(a) 

Figure f: (a) < {jjip > vs. f3 for L = 6** (diamonds) and 
(b) Susceptibility of < ipip > vs. /3 for L = 12^ X 4. 



(b) 



12^ X 4(circlcs), 



4 RESULTS 

Firstly we determine phase transition point /3c at /i = on 12'^ x 4 lattice. 

Chiral condensate < ip-ip > versus (3 is shown in figure 1 (a) . Chiral restora- 
tion occurs at around /3 « 1.4. In order to determine the critical point more 
precisely, susceptibility of the condensate is analyzed as shown in figure 1(b). 
By this analysis, the critical point for m = 0.07 is determined as 



Pc 



1.41 ±0.03 



(5) 



This value is slightly smaller than that of finite size crossover transition 
point f3c = 1.54 on 6* lattice. 

Nextly, we study quark condensate and pseudoscalar meson mass at small 
but finite /i. We stay in confinement phase (/? = 1.1). In figure 2(a), we plot 
the < 'ipip > as a function of /i ( 0.0 ^ 0.2 ) at several m (0.05,0.07,0.1). As 
shown in the figure, the condensate decreases with the chemical potential. But 
it varies very weakly around /i = 0. 

Pseudoscalar meson mass is extracted in a standard way. Typical examples 
of pseudoscalar correlator is shown in figure 2(b). In the present region of 
study, single pole fit gives reasonable fitting and results for pion mass squared 



are shown in figure 3(a). In order to examine the influence of chemical potential 
on hadron mass, we examine pseudoscalar meson consisting of light and heavier 
quarks, figure 3(b) shows their masses at small /z. As shown in the figure, they 
are very stable and no appreciable dependence on chemical potential is seen in 
this region. 
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Figure 2: (a) < 'ip'ij) > vs. ^ for niq = 0.05,0.07,0.1; (3 = 1.1. 

(b) correlation function oi < qq > vs. /3 for /3 = 1.1. fi = 0.1, ruq = 0.05, 0.07, 0.1. 
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Figure 3: (a) M^^ ^g^ „^ foj, p = n^ /.j = ; (b) M^ vs. /i for /3 = 1.1, 



5 SUMMARY AND CONCLUSIONS 

In this paper, we study 8 flavors SU(2) QCD at finite temperature and smafl 
jjL on 12'^ X 4 and 12^ x 24 x 4 lattices. The chiral transition at /i = is 
determined to be (3c = 1.41. Chiral condensate decreases with /i but variation 
is very weak around /i = 0. Mass of pseudoscalar meson consisting of light and 
heavier quarks is stable against chemical potential. 
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